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Abstract. - We show that space-charge instabilities (electric field domains) in semiconductor 
superlattices are the attribute of absolute negative conductance induced by small constant 
and large alternating electric fields. We propose the efficient method for suppression of this 
destructive phenomenon in order to obtain a generation at microwave and THz frequencies in 
devices operating at room temperature. We theoretically proved that an unbiased superlattice 
with a moderate doping subjected to a microwave pump field provides a strong gain at third, 
fifth, seventh, etc. harmonics of the pump frequency in the conditions of suppressed domains. 



Introduction. - There exists a strong demand for miniature, solid-state, room tempera- 
ture operating sources and detectors of THz radiation (0.3-10 THz). The need is caused by a 
rapid progress of THz sciences and technologies ranging from the astronomy to the biosecurity 
PP. Semiconductor superlattices (SLs) operating in the miniband transport regime, are 
interesting electronic devices demonstrating properties of both nonlinear and active media. 
Nonlinearity of voltage-current (UI) characteristic of SL gives rise to a generation of harmon- 
ics of microwave and THz radiation 012]. On the other hand, Bloch oscillations of electrons 
within a miniband of SL cause an appearance of negative differential conductance (NDC) for 
dc fields (voltages) larger than the critical Esaki-Tsu E c |2J. Remarkably, the static NDC 
should be accompanied by a small-signal gain for ac fields of very broad frequency range from 
zero up to several THz |2J. The feasible device, employing such active media properties of 
dc-biased SL at room temperature, is known as THz Bloch oscillator. However, an existence 
of static NDC causes space-charge instabilities and a formation of high-field domains inside 
SL [Hj. These domains are believed to be destructive for the Bloch gain. Recent top class 
experiment demonstrates only decrease in THz absorption, but still not gain, in an array of 
short dc-biased SLs [7j- Therefore, it is important to consider modifications of this canonic 
scheme of generation in order to suppress electric domains but preserve the high-frequency 
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gain in SLs. Recently, the experimental realizations of SL oscillators with ac bias have been 
reported 00. 

In the present Letter, we analyze theoretically the scheme of high-frequency superlattice 
oscillator where instead of dc bias a quasistatic ac electric field is used. Because for typical SLs 
the characteristic scattering time at room temperature is of the order of 100 fs, the quasistatic 
interaction of ac field with miniband electrons is well defined for the microwave frequencies. 
SL is placed into the resonator tuned to the frequency of generation. Here the ac field plays 
a twofold role: It suppresses the space-charge instability inside SL and simultaneously pumps 
an energy for generation and amplification at higher frequency. 

The main idea of our scheme for suppression of space-charge instability is very simple. 
Suppose that the frequency of the pump field lu is larger than the inverse characteristic time of 
domain formation tTt . This requirement, LOTdom S> 1, imposes the lower limit for the values 
of pump frequency lu in our scheme, which, however, still easily belongs to the microwave 
range. Then, we show that the electric instability in the ac-driven superlattice arises if the 
slope of dependence of the time-averaged current on the applied dc bias is negative (NDC for 
the time-average current). We note that the only possibility to have NDC for the time-average 
current in the limit of small dc bias is related to the existence of absolute negative conductance 
(ANC), when the direction of dc current is opposite to its direction in a normal conductor 
(fig. HI). ANC is known in various systems |10| . but this nonequilibrium phenomenon requires 
very special conditions, which as rule are quite difficult to satisfy. In particular, only a few 
combinations of quasistatic ac fields can induce ANC in SL. Thus, a high electric stability 
of the oscillator's scheme is based on a simplicity of ANC state suppression in unbiased or 
weakly dc-biased SLs. 

We prove that within this stabilization scheme the only way to get gain for both small- and 
large-signals, is to tune a resonator to some particular harmonics of to (fig. 0) ■ The physical 
processes in SL contributing the the net gain are the parametric amplification of the field 
in resonator, seeded by the the frequency multiplication of the pump field, as well as the 
nonparametric absorption controlling the suppression of domains. 

Suppression of domains. We consider a response of electrons, belonging to a single 
miniband of SL, to the action of strong pump field E p — Ed c + E u cos cut, where dc bias Ed c is 
small and ac field is quasistatic lot < 1 (r is the intraminiband relaxation time). This pump 
field creates the voltage U p = E p L across SL of the total length L. We suppose also that 
due to a nonlinear character of miniband transport, a signal with amplitude E\ < E^ and 
frequency lu\ > lu is generated. The total voltage across SL U(t) is 



For lot < 1, the dependence of the current through SL I(t) on the ac voltage U(t) can be well 
described by the Esaki-Tsu formula 



where U c = (hL)/ear is the critical voltage (a is the SL period), and I p = I(U C ) is the peak 
current 0. Now we turn to the derivation of electric stability conditions of SL device. For 
lot < 1 a space-time evolution of the electron density p(x,i), the current density j(x,t) and 
the field E(x,t) inside SL, driven by the given voltage U(t), can be well described by the 
drift-diffusion model Set of equations consists of the current equation j = epV(E) — 

D(E)dp/dx, the Poisson equation dE/dx = Ane^ 1 (p — N) with N being the doping density, 
the relation to the applied voltage U(t) — J Q E(x, t)dx, and the continuity equation edp/dt + 



U = U p (t) + U s (t) = (U dc + If,, cosujt) + Ui cos Lu x t. 



(1) 



I = 2I P (U/U C )/[1 + (U/U C )% 



(2) 
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Fig. 1 - [color online] Dependence of time-averaged current (7) on dc voltage Udc for the pump 
Uui/Uc = 9 and for different strengths of signal Ui at frequency 7lo: Positive slope for Ui/U c = 0.9 
(green online), negative slope and ANC for Ui/U c = 6.3 (red online). Esaki-Tsu characteristic 
(U u =0) is shown as dashed line. 

dj/dx = • Dependence of the electron velocity V on the local field E(x, t) is determined 
by the Esaki-Tsu formula as V(E) = 2V p (E/E c )/[\ + (E/E c ) 2 ] with E c = h/ear and V p = 
V(E C ). Following the Einstein relation the diffusion coefficient is D = ksTV(E)/ E. If the 
amplitude of U(t) > U c , space-charge instabilities periodically arise during the part of period 
of ac voltage T = 2-k/uj when the SL is switched to the state with NDC. 

We suppose that condition uJTdom 3> f is satisfied. This inequality guarantees that a charge 
accumulated during T is always limited. In this situation, fluctuations of j and p are small 
and we can use linear analysis investigating stability of the system. The characteristic time 
T dom is of the order of the dielectric relaxation time ~ (w^t)" 1 (uj p i — (47re 2 N/ errio) 1 / 2 is the 
miniband plasma frequency, N is the doping density, tuq is the electron mass at the bottom 
of miniband, e is the averaged dielectric constant of SL) |12| . In more accurate approach 
Tdom depends on the amplitude of ac field. Our numerical simulations demonstrate that for 
typical semiconductor SL with N ~ I0 16 cm -3 accumulated charge is small (wr^om ^ 1) if 
w/27r > 100 GHz. Decreasing doping density N one can satisfy the condition using lower 
pump frequencies. 

Linearizing the drift-diffusion model equations, we find that small fluctuations of space- 
charge with a long wavelength will not grow, if the dependence of the time-averaged current 
(I) on the dc voltage Udc has positive slope, i.e. 

A={I'[U(t)})>0, (3) 

where averaging (...) is performed over period T, prime means first taking the derivative in 
respect to Udc and then assuming Udc — > 0. Here and in what follows, we will consider the 
current / and the voltage amplitudes U u in units of I p and U c , correspondingly. Combining 
eqs. © and © for the quasistatic pump U = U p (t), we get A p = A(U P ) = (I + U 2 )~ 3/2 , 
which is always positive for all U u . Therefore, SL driven by a monochromatic and quasistatic 
ac field is stable against small fluctuations of charge or field. 

We also considered the influence of signal fields on the stability of the system calculating 
numerically the increment A for U(t) given by eq. Q with u)\ = mot. For n — 3 we found that 
A > for all U\ < 10. Nevertheless, the UI characteristics with ANC, and as a consequence 
with A < 0, can exist for some particular values of U u and U% in the cases n = 5 and n = 7 
(fig. [!}. However, for the practically important range U\ < {/*' (see next section), we still 
found no ANC. 

We should distinguish our scheme of electric stabilization based on avoidance of ANC in 
unbiased SLs from the free of domains Limited Space Accumulation (LSA) regime, which is 
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Fig. 2 - Schematic representation of the domainless oscillator: Generation at n = 3, 7 . 



well-known in the Gunn diodes (see J3| f° r the theory and ^2] f° r the theory and experiment) 
and can be applied to SLs as well 15 . This LSA regime works only for a large dc bias and 
for a large ac field in the resonator. In contrast, within our scheme it is possible to avoid a 
formation of domains even for small-amplitude fields in resonator. 

Small-signal gain. - In the oscillator, SL is placed in a resonator providing a feedback 
only for the signal with frequency u>\. Absorption of the signal field in SL is defined as 
A = (I[U(t)] cos(wit)) t , where U(t) is given by and averaging (. . .) t is performed over a 
period that is common for both U p (t) and U s (t). Gain corresponds to A < 0. 

In our case the pump is always quasistatic (wt < 1), but an interaction of the signal 
field with miniband electrons can be both quasistatic for sub-THz frequencies or dynamic 
for oj\t > 1. Accordingly, to find I(U) for lot < 1 and io\t < 1 we can use the Esaki-Tsu 
eq. ©, but have to use the formal exact solution of Boltzmann transport equation ^B] if 
io\T > 1. Surprisingly, in our situation both these approaches give quite similar physical 
results, therefore we focus mainly on more simple quasistatic calculations and only briefly 
discuss a gain for signal with lo\t > 1. 

For a small signal U\ <C U u , we expand I(U) w I(U P ) + I'(U P ) X U s and substitute this 
expression in the definition of absorption A. We see that because of averaging over time the 
absorption A is strongly dependent on the ratio of frequencies u>x and to. First, if to and 
iO\ are incommensurable, we have A — A inc = (I'(U p ))Ui/2. Obviously, for the Esaki-Tsu 
dependence A inc — A p Ui/2 > 0. Moreover, A > also for the probe with lo\t > 1. Second, 
if u>\ = nu>/2 (n is an odd number), we found that always A > for A p > 0. Thus, no 
gain at incommensurate or half-integer frequencies is possible in the conditions of suppressed 
domains. 

Next, we consider the case when the probe frequency is uj\ — nu (fig. 0). In this case, 
the total absorption A is the sum of three terms: Defined above Ai nc , as well as Ah = 
{I (Up) cos(nwi)} and A ca h = (I'(U P ) cos(2nujt))Ui/2, which can be called the incoherent, the 
harmonic and the coherent absorption components, correspondingly |17) . 

We see that while A co h and Ai nc are dependent on both U& and U\ , the term Ah is a func- 
tion of only pump strength U u , and therefore Ah gives the main contribution to the absorption 
of a weak probe. For the Esaki-Tsu dependence we find A h = (-l) fe [2(6 - l) 2k+1 ]/bU% k+1 , 
where b = (1 + U 2 ) 1 / 2 and 2k + 1 = n (k = 1,2,3...). This equation describes just odd 
harmonics of the current (cf. Eq. 17 in 0]); in the limit of weak pump U w <C 1 it takes 
familiar form oc t/™ 0. Importantly, Ah is negative for the odd values of k (fig. OH)- There- 
fore, generated harmonics with n — 3, 7, . . . can provide seeding gain for an amplification of a 
probe field. In quasistatic picture, it has no threshold in the amplitude of pump U u . Thus, 
if the pump amplitude is less than the critical Esaki-Tsu voltage, the gain at 3w will be not 
accompanied by the space-charge instabilities in SL. 

Nonparametric effects in absorption are described by the term Ai nc cx A p . It is always 
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Fig. 3 - [Color online] (a) Amplitudes of the 3d (red online), 5th (blue online) and 7th (green online) 
harmonics of the current, Ah, as functions of the pump amplitude Uu. (b) Incoherent Ai nc (black) 
and absolute values of coherent — A co h components of small-signal absorption at 3u>, 5u and 7u>, as 
functions of pump [7 W . 



positive because of space-charge stabilization. The term A co h describes a parametric ampli- 
fication of the probe field due to a coherent interaction of the pump and the probe fields in 
SL. Numerical calculation of the integral demonstrates that A co h is always negative for odd 
n = 3, 5, 7; its absolute value is larger than Ai nc for a large enough pump strengths (fig. dp)- 
In the limit of weak pump [/„ < 1, we find A co h oc U^ n U\. Now we can describe the ampli- 
fication of a weak signal at 3a; under the action of a weak pump. Third harmonic of a weak 
pump is generated at the cubic nonlinearity of SL characteristic, then it can get a seeding 
gain at the same nonlinearity {Ah oc U^), but the next step of parametric gain, A co h, uses 
already the seventh order (~ U 7 ) nonlinearity in the I(U) dependence. 

We should underline that the sign of generated current at the frequency of harmonic plays 
a very important role in the process of amplification. Really, let us compare the action of 
SL oscillator tuned either to 3ui or to 5w modes (fig. El- Net gain at 3oj is quite large, what 
can provide an effective operation of the oscillator. On the other hand, a radiation at 5th 
harmonic also can be generated in SL. However, in contrast to the case of the 3rd harmonic, it 
cannot be further amplified (because Ah > and Ah > \A co h\) and field in the cavity, tuned 
to 5u>, eventually evolves to zero. 

Large-signal gain. - We calculate the large-signal absorption A(U U , U±) using two meth- 
ods: (i) simple ID calculations with Esaki-Tsu characteristic (eqs. and Q) and (ii) 3D 
single-particle Monte Carlo computations with an account of electron scattering at optical 
and acoustic phonons ^HJ. For the Monte Carlo computations we consider GaAs/AlAs SL of 
the period a — 6.22 nm and the miniband width A = 24.4 meV. Static UI characteristic of 
this SL can be well described by the Esaki-Tsu formula with E c — 4.8 kV/cm and r = 220 fs. 
We took w/2-k = 100 GHz. These are our default parameters. 

The dependence A{Uu) for 3w-generation and for different values of the relative probe 
amplitude, 77 = Ui/U u , are shown in fig. 01 Results of simple ID theory and 3D Monte Carlo 
simulations are in a good agreement. For rj = and for small 77 <C 1, the dependence of A^U^) 
follows to the corresponding dependencies for the seeding gain, Ah(U w ), and the small-signal 
gain. With a further increase of ?/, the gain decreases and finally the absorption becomes 
zero for some 77 = j]q: A(U u ,r]o) = 0. The value t]q = [llf/U^] 2 determines the maximal 
device efficiency for the given pump (here U[ ' is the voltage amplitude corresponding to the 
stationary amplitude of ac field in the resonator). Inset in fig. E| shows the smooth dependence 
of 770 on the pump amplitude. Generation with low efficiency (less than 5%) is possible even 
without static NDC, for U u < 1. The efficiency reaches its maximal value of 23% at U u ~ 4. 
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Fig. 4 - The dependence of absorption at 3ui on the pump amplitude Uu, for different relative probes, 
n. Inset: Maximal relative amplitude of the signal allowing gain, 770, as a function of the pump. 
Calculation with Esaki-Tsu characteristic (solid) and Monte Carlo technique (symbols). 

Fig. 5 - Regime of suppressed domains in SL oscillator. Relative reduction of gain at 3d harmonic 
S vs voltage amplitude U u for uj/uid = 0.1 and 77 = U\/U u = 0.1. Inset: S as a function of the ratio 
uj/uJd for Uuj/Uc = 2, rj = (squares) and r\ = 0.4 (triangles). 

For 7w-oscillations the maximum of t]q is 4%. 

Numerical evidence. - We directly calculate the influence of inhomogeneous distributions 
of space-charge on the gain in SL. We numerically solve the drift-diffusion model for the 
time-dependent voltage U(t) (eq. QJ) with Ud c = and determine the absorption at the nth 
harmonic as Ad = (((J(x, i)/jo) cosnwi}), where jo = eV p N and averaging ((. . .)) is performed 
both over the period T and the length L. In the computations a spatial fluctuation of charge 
distribution was caused by a local Gaussian-profiled deviation of E c (or r). We observed a 
periodic spatial accumulation of charges during a part of period T. We computed the relative 
decrease in gain, 8 = (A — Ad) /A, for different harmonics n and pumps t/ w . Dependence 
of 8 on U u , for SL of 130 periods and for uj/uJd = 0.1 (uid = v^t), is shown in fig. |SJ In 
fig. El the reduction of gain due to charge accumulation is 2% for optimal = 4 and it is 
less than 8% overall. For our default parameters, the value uj/uid = 0.1 corresponds to the 
doping N = 2 x 10 16 cm~ 3 . In this situation uJTdom — L We found that with an increase 
of uj/uid oc UJTdom j the value of 8 quickly decreases (see Inset in Fig. |SJ|. Therefore, for 
higher frequencies or lower doping (uid oc N) providing uijuid > 0.1, the influence of charge 
accumulation on the gain is practically negligible. 

Additionally, we also have performed the ensemble Monte Carlo simulations to modell 
the oscillator with account of space-charge dynamics in SL. In general, these 3D simulations 
do confirm our main conclusions on the suppression of domains made in the framework of ID 
drift-diffusion model. 

Conclusion. In summary, we derived the necessary and sufficient conditions for the 
suppression of space-charge instabilities in the high-frequency SL oscillators with a microwave 
pump. Our findings can provide simple tests to distinguish domainless and domain-mediated 
regimes of generation in these oscillators, such as a measurement of the dependence of the 
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time-average current on the dc voltage or a comparison of operation at 3lu and 5w. 

The domainless oscillator can operate at THz frequencies. Probably the most easily re- 
alizable case is to use a pump with u> near 100 GHz to obtain an output in the important 
frequency band of hundreds of GHz. In order to suppress electrical domains and to use res- 
onator with a reasonable Q (e.g., Q < 100 for generation at 3w), the doping of SL should be 
~ 10 16 cm" 3 . 

In the recent experiments [HIE]) a generation at the third harmonic of a strong microwave 
pump of 100 GHz have been reported in SL device. This generation has been first attributed 
jS] to the process of periodic formation of electric domains, but later it has been re-attributed 
to a parametric generation in [5]. In (8JE] a heavy doped (N — 10 18 cm' 3 ) SL was used, and 
therefore this device certainly operated in the regime with uT dom <C 1, which is opposite to the 
condition of limited charge accumulation. In this situation, a sufficient charge accumulation 
during every cycle of pump field arises, and therefore SL can switch to ANC state, as has 
been shown for the case of microwave-driven Gunn diode in [20| . We also saw such kind of 
ANC states in our simulations of SL oscillator within the drift-diffusion model. This ANC 
induced by the accumulated charge completely alters electric stability of SL device and results 
in an operational mode that is quite different from the mode based on an avoidance of ANC. 
Detailed theoretical analysis of the experiments [HIE] goes beyond the scope of present paper. 

* * * 

This research was supported by Academy of Finland (grants 202697 and 109758), Emil 
Aaltonen Foundation, Programme of RFAE (grant 4629) and AQDJJ Programme of ESF. 

Erratum. - Additional modelling of a high-frequency generation with consideration of a 
resonator tuned to nu> showed that the phase difference 4> n between the pump field (lo) and 
the signal field in resonator (nu) can be either for n — 3,7 or -k for n = 5,9 within the 
quasistatic approximation. Condition Ah > determines the phase shift (j)^ — n between 
the pump field (to) and the signal field in resonator (5cj), but it does not prevent a growth 
of this mode, in contrast to our previous statement. Therefore, the amplitude of field E\ 
in an ideal resonator tuned to 5w in fact can grow until it reaches a stationary value 
corresponding to zero value of the total absorption A(E LU , Ef) = 0, which is now defined 
as A = (I[E(t)]cos(nujt + (j) n ))t for the total field E(t) = E^cos(ujt) + EiCOs(nuit + 4> n ). 
The maximal efficiency ?7q = [Ef / EJ\ 2 is 10%. Next, following the necessary condition for 
suppression of domains [eq. (3)] the generation at 5w is still possible in domainless mode. 
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